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'~{ _ Abstract 

\ In this paper, a backward Euler method is discussed for the equations of motion arising in the 2D Oldroyd 

model of viscoelastic fluids of order one with the forcing term independent of time or in L°° in time. It is shown 
■ that the estimates of the discrete solution in Dirichlet norm is bounded uniformly in time. Optimal a priori 

error estimate in L^-norm is derived for the discrete problem with non-smooth initial data. This estimate is 
shown to be uniform in time, under the assumption of uniqueness condition. 
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1 Introduction 

In this paper, we consider fully-discrete approximations to the equations of motion arising in the Oldroyd fluids 
^ ■ (see Oldroyd 15 , Oskolkov[16 ) of order one: 

^! (1.1) ^ + u- Vu-/.tAu- / P{t-T)Au{T)dT + \7p = {. mn,t>0 

OO ■ with incompressibility condition 

o 

CN ■ (1.2) V • u = 0, on f], i > 0, 



and initial and boundary conditions 

(1.3) u(x,0)=uo in 17, u = 0, on Srj, i > 0. 

Here, $7 is a bounded domain in with boundary /i — 2kA~^ > and the kernel /3(t) = 7exp(— where 
7 = 2X~'^{v — k\~^) > and 5 = A^^ > 0. Further, f and Up are given functions in their respective domain of 
definition. For more details, we refer to pQ and [IS]. 

There is considerable amount of literature devoted to Oldroyd model by Oskolkov, Kotsiolis, Karzeeva, Sobolevskii 
etc, see [U m [m [m [16] and recently by Lin et al. [9l[10l[24], Pani et al. [19l|20], Wang et al. [25], and references, 
therein. A detailed report on the continuous and semi-discrete cases can be found in 

Literature for the fully-discrete approximations to the problem (|l.l|) - ()1.3p is, however, limited. In 2 , Akhmatov 
and Oskolkov have discussed stable and convergent finite difference schemes for the problem (|l.ip - (|1.3|) . Recently 
in [20) . a linearized backward Euler method is used to discretize in temporal direction and semi-group theoretic 
approach is then employed to establish a priori error estimates. The following error bounds are proved in |20] for 

tn>0 

l|u(i„)-U"|| <Ce-"*"fc 

and 

||u(t„) - U"||i < Ce-"*"fc(i-i/2 + log i) 
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for smooth initial data and for zero forcing term. Here, k is the time step and U" is the finite difference approx- 
imation to u(t„), when modified backward Euler method is applied in the temporal direction. Recently Wang et 
al. [5S] have again applied backward Euler method for the problem (|l.ip - (|1.3p . with smooth initial data, when 
the forcing function is non-zero. They have used energy arguments along with uniqueness condition to obtain the 
following uniform error estimates: 

||u(i„) -U"|| < C{h'^ + k) 

and 

ri/2||u(t„)-U"||i <C(/i-)-A:), 

where T(t„) = min{l,t,i} and h is the mesh size, again with smooth initial data. 

Our present investigation is a continuation of [S], where a priori estimates and regularity results have been 
established, which are uniform in time under realistically assumed regularity on the exact solution and when 
f , ft G L°°(L^). Error estimates for semi-discrete Galerkin approximations have been shown to be optimal in 
L°°(L^)-norm for non-smooth initial data. Further, uniform (in time) error estimates under uniqueness condition 
are also established. 

In the present article, we discuss backward Euler method to discretize in the temporal variable and Galerkin 
approximations to discretize spatial variables for approximating solutions of the problem (ll.ip - (|1.3p . Our main 
aim, in this work, is to present optimal error estimate for the backward Euler method, when the initial data is 
non-smooth, that is, Uo G Ji. The main results of this paper are follows: 

(i) Proving uniform bound in time in the Dirichlet norm for the solution of the completely discrete backward 
Euler method. 

(ii) Deriving new estimates which are valid uniformly in time for the error associated with discrete linearized 
problem 

(iii) Establishing estimates for the error related to nonlinear part in which the error constant depends exponentially 
in time and thereby, making final error estimate in the velocity to depend on exponentially in time. 

(iv) Proving optimal error estimates for the velocity in L^-norm which arc uniform in time under the uniqueness 
assumption. 

To prove estimate in the Dirichlet norm for the discrete solution which is valid for all time, the usual tool, in the 
case of the Navicr-Stokes equations, is to apply discrete version of uniform Gronwall's Lemma. Now for proving (i), 
it is difficult to apply uniform Gronwall's Lemma due to presence of the discrete version of integral term. Therefore, 
a new way of looking at the proof has helped to achieve (i), see; Lemma 4.3. For (ii) — (Hi), we observe that there 
are difficulties due to the non-linear term along with the presence of integral term in the case of non-smooth initial 
data. For example, the preliminary result (L°°(L^) estimate) is sub-optimal due to non-smooth initial data (see; 
Lemma l5.2l) . In order to compensate the loss in the order of convergence, a more appropriate tool is to multiply by 
t. But, unfortunately, it fails here due to the presence of the integral term (or the summation term). To overcome 
this difficulty, we modify some tools from the error analysis of linear parabolic intcgro-differential equations with 
non-smooth data (see; [171 IHl [23]) to fit into the present problem and also a special care is taken to bound the 
nonlinear term. Our analysis makes use of the solution, say; V" of a linearized discrete problem (see; (5.5)) as an 
intermediate solution. Then, with its help, we split the error: uJJ — U" at time level t = tn, where uJJ — u/i(i„) is 
the solution of the semi-discrete scheme aX t = tn and U" is the solution of the backward Euler method, into two 
error components: one in ^" := — V", which denotes the contribution due to the linearized part (see; (15. 6p ). 
and the other in ry" := U" — V", which is due to the non-linearity (see; (|5.6p ). Using a backward discrete linear 
problem and duality type argument along with an estimate of where 

n 

an L^-estimate of ^" which is valid for all time is derived, refer to Theorm 5.1. For estimate of rj", we employ 
negative norm estimate and estimate of r)" and obtain estimate which depends on exponentially in time, see; 
Lemma 5.9. Thus, one of the main result for nonsmooth initial data that we have derived in Theorem 5.2 is as 
follows: 

(1.4) Mtn) - V^ < KTtn''\e + k{l + \og\f'^), 
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where Kt depends exponentially on time. Finally for the proof of (iw), a careful use of the uniqueness condition, 
it is also shown that the error estimate (jl.4|) is valid for all time. 

The remaining part of this paper is organized as follows. In Section 2, we discuss some notations, basic 
assumptions and weak formulations. In Section 3, a semidiscrete Galerkin method is discussed briefly. Section 4 
is devoted to backward Eulcr method. Optimal and uniform error bounds are obtained for the velocity when the 
initial data are in Ji. Finally, we summarize our results in the Section 5. 



2 Preliminaries 

For our subsequent use, we denote by bold face letters the R^-valued function space such as 

= [Hl{Vt)f, = [L^{n)f and H™ = [H"\9.)f , 
where H"'-{n) is the standard Hilbert Sobolev space of order m. Note that Hj is equipped with a norm 

Further, we introduce some more function spaces for our future use: — 

Ji = {0 e Hj : V • ^ = 0} 
J = {0 £ : V • = in O, • n|aj2 = holds weakly}, 

where n is the outward normal to the boundary 917 and cf> ■ n\gfi = should be understood in the sense of trace in 
H~i/2(9f2), see [22] • Let H"^/R be the quotient space consisting of equivalence classes of elements of ff™ differing 
by constants, which is equipped with norm = ||p + c||,„. For any Banach space X, let LP(0,T; X) denote 

the space of measurable X -valued functions 4> on (0,r) such that 



[ \\cl){t) fx dt < oo ifl<p<cx), 
Jo 



and for p = oo 

ess sup ||c/)(i)||x < cxo if p ^ og. 



0<t<T 

Through out this paper, we make the following assumptions: 

(Al). For g € L^, let the unique pair of solutions {v G Ji, 9 G L^/R} for the steady state Stokes problem 

-Av + Vg = g, 
V-v = in n, v|an = 0, 

satisfy the following regularity result 

||v||2 + ||<z||ffV«<^l|g|l- 

(A2). The initial velocity Uq and the external force f satisfy for positive constant AIq, and for T with < T < 00 
uo G Ji, f,ft G L°°(0,r;L2) with ||uo||i<Afo, sup {||f ||, ||ft||} < Mq. 

0<t<T 

For our subsequent analysis, we present the following Lemma, which can be seen as a discrete version of Lemma 
2.2 from 19 . 

Lemma 2.1. Let g^, 0' G M, 1 < i < ri, n G N and < k < 1. Then the following estimate holds 

i—l J — 1 i—1 i—1 
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3 Semidiscrete Galerkin Approximations 

From now on, we denote h with < ft, < 1 to be a real positive discretization parameter tending to zero. Let Hh 
and Lfi, < h < 1 he two family of finite dimensional subspaces of Hj and L^, respectively, approximating velocity 
vector and the pressure. Assume that the following approximation properties are satisfied for the spaces and 
Lh-. 

(Bl) For each w G Hq n and q E / R there exist approximations ihw g H/j and jhq G Lh such that 

||w-iftw|| + ft||V(w - iftw)|| < ifo/i^||w||2, Ik- j/iglUvfl < KQh\\q\\Hi/B,- 
Further, suppose that the following inverse hypothesis holds for w/j G H^: 

(3.1) ||Vw,,|| < A'o/i-i||w,,||. 
For defining the Galerkin approximations, set for v, w, G Hj, 

a(v,0) = (Vv,V</.) 

and ^ ^ 

6(v,w, 0) = -(v Vv^r, 0) - -(v- V0,w). 

Note that the operator &(•, •, •) preserves the antisymmetric property of the original nonlinear term, that is, 

K^h, Wft, Wft) =0 Vv,i, v^f/, G H/i. 

Now, the semidiscrete Galerkin formulation reads as: Find u.h{t) £ H/j and Ph{t) G Lh such that u;i(0) ~ Uq/i and 
for t > 

{uhu<t}h)+ t^O'{nh,(t)h) + &(u/!.:Uft>h) +a(u/,,/3,0;,) - (p/,, V • 0^) = (f,0^), 

(3.2) (V-u;„X/.) =0, 

for 0^ G Hft,, x/t ^ Lh- Here Ug/i G H/^ is a suitable approximation of Ug G Ji and 

(3.3) MhAt) = / P{t~s)vih{s) ds. 

Jo 

In order to consider a discrete space analogous to Ji, we impose the discrete incompressibility condition on H/j 
and call it as Jh- Thus, we define J/,,, as 

J/i = {vh e H,i : (xh, V -Vh) ^0 Vxh e L/i}. 

Note that 3h is not a subspace of Ji. With J/,, as above, we now introduce an equivalent Galerkin formulation as: 
Find Uh{t) G Jh such that u/j(0) = uoh and for t > 

(3.4) (u,it, 0/i) + A'a(u,i, 0;^) + a{uh,i3, <t>h) = -H'^h, u^, 0/J + (f, 0/i) y4>h £ J/i- 

Since is finite dimensional, the problem (j3.4p leads to a system of nonlinear integro-differential equations. 
For global existence of a solution pair of ()3.4p , we refer to [19l . Uniqueness (of p) is obtained on the quotient space 
Lh/Nh, where 

Nh ^ {qn e Lh : (g/., V . 0;,) = 0, V0,, G H,J. 

The norm on Lh/Nh is given by 

II^IIlVJV,. = inf hh + XhW- 

For continuous dependence of the discrete pressure Ph{t) G Lh/Nh on the discrete velocity Uh{t) G J/i, we assume 
the following discrete inf-sup (LBB) condition for the finite dimensional spaces H/j and Lh'- 

(B2') For every qh Lh, there exists a non-trivial function (f>h G H/i and a positive constant Kq, independent of 
h, such that 

|(<7;„V.0J| >i^ol|V0,J|||(Z/.l|LViV.• 
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Moreover, we also assume that the fohowing approximation property holds true for Jh- 
(B2) For every vir G Ji n H^, there exists an approximation r/iW e Jh such that 

||w - r^wll + h\\Viw - rhw)\\ < K^h^Wwy. 

This is a less restrictive condition than (B2') and it has been used to derive the following properties of the 
projection Ph : Jh- We now state without proof these results. For a proof, see ^IJ^. For <p ^ 3h, note that 

(3.5) 110 - P„0|| + h\\VPh4>\\ < Ch\\V4>l 
and for e Ji n H^, 

(3.6) 110 - P„0|| + /i|lV(0 - P„0)|| < Ch^\\Act>\\. 

We now define the discrete operator Ah : H/i H> H/i through the bilinear form a(-, •) as 

(3.7) a(v/i, (pi^) = {-AhVh; 0) yvh,4>h e 

Set the discrete analogue of the Stokes operator A = PA as A;i = PhAh- Using Sobolev imbedding and Sobolev 
inequality, it is easy to prove the following Lemma 

Lemma 3.1. Suppose conditions fBlj and (B2) are satisfied. Then there exists a positive constant K such 

that for V, w, G Hh, the following holds: 



(3.; 



(v • Vw,0)| < K < 



l/2||Vv||l/2||Vw||l/2||A„w||l/2||0||, 

i/2||A,v||i/2||Vw||||0||, 
i/2||Vv||V2||Vw||||0||V2||v0||i/2, 

||Vw||||0||l/2||A^^||l/2, 
||Vw||l/2||A,^w||l/2||^||l/2||V0||l/2. 



Examples of subspaces and Lh satisfying assumptions (Bl), (B2'), and (B2) can be found in [6l|4l[3]. 

We present below, a Lemma, that deals with higher order estimates of Uh, which will be useful in the error analysis 

of backward Euler method for non-smooth data. 

Lemma 3.2. Suppose conditions (Al), (Bl), (B2) and (BA) are satisfied. Moreover, let Uh{0) G 3h andt satisfy 
the assumption (A.3). Then, Uh, the solutions of the semidiscrete Oldroyd problem Jg.^p satisfies the following a 
priori estimates: 



(3.9) 
(3.10) 

(3.11) 



T*\\uh\\l + {T*r'-'\\nht\\l<K, rG{0,l}, 
e^''%T*Y{s)\\uhs\\lds<K, rG{0,l,2}, 







-2at 



^T*Y+\s)\\uhss\\l-i ds<K, rG {-1,0,1}, 



where (T*){t) = niin{l,i}, a{t) — T*(t)e and K depends on the given data, but not on time T. 

Proof. The estimates (j3.9|) - p.l0p can be proved as in the continuous case, see [8]. For the final estimate, we 
differentiate (|3.4p to find that, for 0,^ G J/i, 



(u/itt,0^) +/ia(u;»t,0/i) + /3(O)a(u/i,0,J - (5 / /3(i - s)a(u,,(s), 0,J ds 

(3.12) = -fe(u/it,Uh,0,,) -fo(u,i,Uftt,0^) + (ft,0^). 

Taking 0,, = {T*f{t)e'^°'^\ihtt in ((37T2)) . we obtain 

{r*f{t)e'-'Wur + f |((r*)^(t)e^"*||u;,||?) < (a(r*)2(t) + r*(t))e2-||u;,f 



+ 7(r*)^(t)e^"*||u,||2||u,„||+<5(r*)2(t)e2"* 



(3{t - s)\\uh{s)\\2\\uhtt\\ ds 



(3.13) 



+ iT*f{t)e'^''Wb{uht,Uh,Uhu)\ + \b{uh,Uht,Uhtt)\ + ||/t|l llu^ttl 
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Use dS^) to find that 

\b{^ht,Uh,,Uhtt)\ + \b{uh,Uht,Uhu)\ < ^llu/ittf +i4:||u,it||f|ju?i|||. 
Now, using p.9p - p.l0p . we can easily deduce from p. 131) that 

(3.14) {T*)'\\nht\\l + /xe-2«* f\T*f{s)e'"'\\uhss\\' ds < K. 

Jo 

We set = -r*(i)e2"*A-iu,,it in ([5T^ . From ([3^, we see that 



t 

2q«IIY7,, I|2 I ollp ||2 I rU., X\/ I f^lU „^„at|IA ,, ('„MI ^„\2 



and therefore 



+C(/x,7)r*(i)e'"*||Vu„r + 2||ft||^ + C{fi,S){ / /3(i - s)e"*||A;,u;,(s)|| ds) 

Jo 

+Ci^i)T*it)e''-'(\\Vuh\\^Uht\\^ + ||Vu,tf (1 + ||u„||||Vu, 
Integrate with respect to time and multiply by e~^"* to conclude 

(3.15) T*(i)||u,tf + ^e-2"* f T*{s)e^''''\\nhss\tids<K. 

Jo 

Finally, we set — —e^^^Aj^Uhtt in p.l^P and proceed as above to arrive at 

(3.16) ||uw|l% +^^6-2"* [\^"^\\u^,,\\l^ds < K. 

Jo 

This completes the rest of the proof. □ 
The following semi-discrete error estimates are proved in 

Theorem 3.1. Let be a convex polygon and let the conditions (A1)-(A2) and (B1)-(B2) be satisfied. Further, 
let the discrete initial velocity Uoh G J/i with Uq/i — P/iUg, where Uq € Ji. Then, there exists a positive constant C 
such that for < T < oo with t € (0, T] 

||(u - u^,){t)\\ + h\\Viu - u^,){t)\\ < Ce^'hh-'/^ 

Moreover, under the assumption of the uniqueness condition, that is, 

(3.17) ^l|f||oo<l andN^ sup ^^'^'l^'^l 

u,v,weHi(n) l|Vu||||Vv||||Vw|| 

where v = pL + ^ and ||f ||oo ||f ||L°°(o,oo:L2(n)) then we have the following uniform estimate: 

\\{n^nH){t)\\<Chh-^/\ 

4 Backward Euler Method 

For time discretization, we state below some notations. Let fc, < fc < 1, be the time step and let i„ — nk, n > 0. 
We define for a sequence {(/>"},i>o CJ/i, 
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For continuous function v{t), we set v„ = v(t„). Since backward Euler method is of first order in time, we choose 
the right rectangle rule to approximate the integral term in p.4|) as: 

(4.1) q'^icP) = fc V/3„-j0^ « / /3(i„ - s)<l>{s) ds 

3=1 ■'° 

where l3n~j — (3{tn — tj). With Wnj = kp{tn — tj), it is observed that the the right rectangle rule is positive in the 
sense that 



(4.2) k^^lli^W ^kJ2kY.^v<l^<l^'>0, 

1=1 i=l j=0 



For positivity of the rectangle rule with Uno — 0, we refer to McLean and Thomee [2]. Note that the error incurred 
due to right rectangle rule in approximating the integral term is 

(4.3) £:?(</)) / /3(t„ - s)0(s) ds - fc V 

•^0 ,,=1 



n „f - 



< 



^(/3(i„-s)0(5)) 



ds. 



We present here a discrete version of integration by parts. For sequences {ai} and {bi} of real numbers, the following 
summation by parts holds 

i i — 1 

(4.4) k^^ajbj = aSi - ky^{dtaj+i)bj, 
where bi := fc^*^-^ bj. 

We describe below the backward Euler scheme for the semidiscrete Oldroyd problem p.2p : Find {U"}„>o G Hh 
and {P"}ri>i G Lfi as solutions of the recursive nonlinear algebraic equations (n > 1) : 

(9iU", 0,J + A^a(U«, 0J + a(g»(U), 0J = (P", V • 0,J 

(4.5) + (f",0,J-&(U",U",0,J eH,, 

(V-U",xh) = Vx/. ei/., n>0. 

We choose U" — Uoh — P/jUq. Now, for 0^ G J/j, we seek {U"}„>o G Jh such that 

(4.6) (atU", 0J + A*«(U", 0,) + a(<z,"(U), 0,) = (f", 0J - fo(U", U", 0J V./., £ J;,. 

Using variant of Brouwer fixed point theorem and standard uniqueness arguments, it is easy to show that the 
discrete problem (j4.6p is well-posed. For a proof, we refer to [7]. Below we prove a priori bounds for the discrete 
solutions {U"}„>o. 

Lemma 4.1. Let < a < min{5, and ko > be such that for < fc < fco 

l+(^)fc>e"'=. 

Further, let JJ^ — Uq/i = Ph'^o with Uq G Ji. Then, the discrete solution U^, N > 1 of satisfies the following 
estimates: 

N 

(4.7) ||U^||2 + rie-"*"fc^e"*"||VU"f < c(e-"*" ||U0||2 + ||f ||^ 
where ||f||oo = ||f ||l°°(L2) , and 



r,^(e--M-2(^-^)A-) 
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Proof. Setting U" = e"*"!!", we rewrite (|4.6p . for 0,^ G J/j, as 

(4.8) e"*"(atU", 0J + fiaiij", + e-"*"&,(U", U", 0J + e"*"a(g,"(U), 0,J = (f", 
Note that ^ 

On substituting this in (|4.8p and then muhiplying the resulting equation by e^"*^, we obtain 

{dtV\ cf,^) - ( ) (U", 0,J + e-"V«(U", 0J + e-"*"+^&(U", U", J 

n 

(4.9) +7e-"'^' ^ e-('^-")(*"-*')a(U^ 0^) = e~"'=(f", 0^). 

i=l 

Put 0^ = U" in and observe that 



and that the nonUnear term vanishes. Also use ||U"|| < x;-||VU"l| to obtain 

1 / 1 _ p-ak , 

-5,||U"|l2 + (e-"V - ( l|VU"f 

71 

(4.10) +7e-"'''fc^e-(*-")(*"-*')a(U\U") < e-"'^'||f"|| ||U"||. 



i=l 



The right-hand side of (|4.10p can be estimated as 



so as to obtain from (|4.10p 



,1 - e-"''. 



at||u"f +(e-"V - 2{- — r — UrM livu""2 



k 



n 

(4.11) +27e-"'=fcy e-('^-")(*"~*')a(U'',U") < — e-"'=||f"f . 
With < a < min{(5, ^^}, we choose fco > such that for < k < kg 

This guarantees that e~"'^ ^~2{^^—^ — )^r^ — 0- Multiply (|4.1ip by k and then sum over n = 1 to A^. The resulting 
double sum is non-negative and hence, we obtain 

N |jn|,2 N 

(4.12) ||U^f + Tifc ^ II VU"||2 < ||U"||2 + lLl2Sie-"'=A: ^ e^"'". 

n=l ^-^1 n=l 

Note that using geometric series, we find that 

(4.13) k e'"*" = e'^'^^^jf^e^"*" = e2a(fc-fc*)g2at„^ 

n— 1 

for some k* in (0, k). On substituting (|4.13p in (|4.12p . multiply through out by e~"*" to complete the rest of the 
proof. □ 



In order to obtain uniform (in time) estimate for the discrete solution U" in Dirichlet norm, we introduce the 
foUowing notation: 

n 

(4.14) U^ = fc^/3„,U^n>0; U^^ = 0, 



and rewrite (|4.6p . for 0^ G J/i, as 

(4.15) (5tU", 0,J + /ia(U", 0J + 6,(U", U", 0J + a(U^, 0J = (f", </,,). 
Note that 

(4.16) U;^' = fc7U" + e-^^\5}-^, 
and therefore 

(4.17) a^u^^ = i(u;^ - u;-i) = ^(fc7U" + e-^'^u;-! - u;^'!) 

Lemma 4.2. Lei < a < min((5, /^Ai/2), U° = P^uq and ko > be such that for < k < ko 

l+(^)fc>e-. 

Then, the discrete solution U", n > I of |^.^ satisfies the following uniform estimates: 

(4.18) ||u"||2 + -_||VU^f <e-2"*"||U0|p+ -^j ||f||L - Mf^, 



and 



(4.19) fc ^ (A*l|VU"f + -IIVU^IP) < Ml, + — ||f||L = MUl), 

n—m 

where UJ^ is given by ^.14\ l. 

Proof Take 0^ = U" in and from (HTT]) . we find that 

e-^\,,,„ , (l-e-^'^) 



a(U^,U") - ^a(U^,5tU^) + ^— '-\\VVy/^ 

7 K7 



Using mean value theorem, we observe that 



k 

Therefore, we obtain from (I4.15P 



^ '-^Se-^^' >(5e~*^ fc*G(0,fc). 



C i,„^^^iin\ II II 9 ZjUC II II 9 , J- 



(4.20) 5t(||U"f + llvU^f) +^||VU"||2 + -1^ llVU^f < ^||f"f . 

7 7 ^Ai 

As < a < min{(5, /iAi/2}, we now find that 

-<5fe 1 

(4.21) 5((||U"||2 + llVmf ) + 2a(||U"f + ||Vmf ) < — ||f"||'. 

7 '"^ 7 M'^i 

Multiply the inequality (|4.2ip by e"°*"-i for some ao > and note that 

f pOiok _ 1 , 

a,(e"°*"0") = e^"'-^[dtr+ ^ r] 

(4.22) < e^o'-ijat^" + 2a0"}. 



With the assumption on the time step k, that is, < fc < fco, and for given a, we can always choose ao such that 
(4.23) l + 2afc>e"«'=. 

Observe that ao < 2a. Therefore, we obtain from (|4.2ip 



at(e"°*"(i|u"f + ^ — l|vu; 



n 1 1 2 1 



< 



ii pii / / — \ II IIOO 

7 // /iAi 

Multiply by k and sum over 1 to n and then multiply the resulting inequality by e~"°'". Observe that = by 
definition. This results in the first estimate (|4.18|) . For the second estimate (|4.19p . we multiply ()4.20p by fc, sum 
over TO to TO + Z with to, / G AA and use (|4.18p to complete the rest of the proof. □ 

Lemma 4.3. Under the assumptions of Lemma \4-^ the discrete solution U", n > I of J^.g] ) satisfies the following 
uniform estimates: 

— 5k 

(4.24) ||VU"||2 + ^ ||AhU'i||2 < if. 

7 

Proof. Set </),j = — A;iU" in (|4.15p and as in the Lemma l4?2l we now obtain 

5t(||VU"f + ^||A;,U^f) +^||A;,U"||2 + -||VU^f < ||r||||A,U"|| 

(4.25) +|6;,(U",U",-A;,U")|. 
Use Lemma [5TT] to arrive at 

dt (||VU"f + f^llA.U^f ) +|l||A,U"f + ^||A,U^f 

(4.26) <^||f||^ + (^)3M2^||VU"||^ 



For some ao > 0, we find that 



n||2 ^ '^IIA TTn||2 , ^ „2||tt"I|2 



(4.27) ao||VU"r < f l|A.U"||^ + —atW 
Define 

(4.28) 5" = min{ao - (|-)=^A/i\ || VU" f , 2^}. 



With := ||VU"f + £^||A,,U"f , we rewrite as 



(4.29) < -||f ||^ + ^al\\\J^^f = K,,. 

Let {riijigN and {fiijieN be two subsequences of natural numbers such that 

5"' = ao+fiXi - (:^)3M2j|VU"i^ = 26, Vz. 
2/x 



If for some n, 



= ao +A*Ai - (A)3Afi\||VU"f = 25 



then without loss of generality, we assume that n £ {rii} so as to make the two subsequence {rii} and {rii} disjoint. 
Now for m,l e N, we write 

m+l ™ii ™i2 

n=m n—rni n—mi 

(4.30) = fc 5] iao+f^Xi - i-fM!,\\yiJ-\n+k ^ 25. 
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Here, mi, m2, • • • , mjj G {rii} n {m, m + 1, • • • , m + Z} and mi, 7712, • • • , mz2 £ {rii} n {m, m + 1, • • • , m + such 
that ^1+^2 — I + 1. Note that li or Z2 could be 0. Using Lemma we observe that 



m + l ^ ^9 + / 



Q """ ' Q3A//^2 Q3A,f2 

^ n—m n—m 

Therefore, from ()4.30p . we find that 

m-\-l 
n—m 

We choose ao such that (A;^i)(ao 4- — i^i2(^i) ~ 2S{kli) to arrive at 

m-\-l 

(4.31) kJ29"> 25ii+i. 

n—rn 

By definition of 5", we have equality in (|4.31l) and in fact, 5" — 26. Now from (|4.29p . we obtain 

dtE" +2dE'' <Kii. 

As in ()4.22p . we can choose < aoi < a < ^ such that 

5^(gaoit„^n) < e"oi*"-i(at£;" + 2(5£;") < i^iie""i*"-i. 

Multiply by fc and sum over 1 to n. Observe that E'^ = ||VU"|p. Finally, multiply the resulting inequality by 
g-«oi«„ to find that 

< e~"oi*"||vu°f + 

This completes the rest of the proof. □ 
Remark 4.1. As a consequence of the Lemma 4.3, the following a priori bound is valid: 

(4.32) T*(t„)||A;,U"||2 <i^. 

5 A Priori Error Estimate 

In this section, we discuss error estimate of the backward Eulcr method for the Oldroyd model (|l.ip - (|1.3p . For the 
error analysis, we set, for fixed neN, l<n<A^, e„= U" — Uh{tn) = U" — ujj. We now rewrite p.4p at t = t„ 
and subtract the resulting one from (|4.6I) to obtain 

(5.1) (dten, (t)h) + Ma(e«, 0J + a{q';{e), 0J = ^"(u,,)(</.J + 4'(u,O(0J + K'fXcjy^X 

where, 

£;"(u,O(0,J = «„<AJ-(9t<,0J = «„<AJ-i^*" (u,„0Jds 

1 



(5.2) " 2fc7, (t-i«-i)(u,.tt,0Jdi, 

(5.3) e"(u,,)(0^) = a{uh,fi{tn),<t}h)ds - a{q';:{\ih),ct)h) = a{£'i^{^h),(l)h), 
and 

A;i(0J = &«, <, 0J - 6(U", U", 

(5.4) = -6«,e„,0,J - 6(e„,<,(/),J - 6(e„, e„, 0^). 

In order to dissociate the effect of nonlinearity, we first linearized the discrete problem (|4.6p . and introduce 
{V"}„>i G as solutions of the following linearized problem: 

(5.5) (a^V", 0J + A^a(V", 0,) + a(q,"(V), 0J = (f", 0J - ^(uJJ, <, 0,) £ J,., 
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given {U"'}n>i S Jh as solution of ()4.6p . It is easy to check the existence and uniqueness of {V"}„>i €E J/i. 
We now spht the error as: 

(5.6) e„ := U" - < = (U" - V") - - V") =: ry,, - 
The following equations are satisfied by and r/^, respectively: 

(5.7) (at^„, <AJ+Aia(^„, <AJ + 0h) = -i?"K)(0J - £^(u„)(c/,,) 
and 

(5.8) (5tJ7„, ^,J+Ma(»7„, </>/J + a(9"(^), 4>h) = Ki4>h)- 
Below, we prove the following Lemma for our subsequent use. 

Lemma 5.1. Let r, s G {0,1}, = min{l,ti} and a as defined in Lemma 4-L Then, with E"' and defined, 
respectively, as I15.2\) and 115. 3\) . the following estimate holds for n = I, ■ ■ ■ , N and for {(p\}i in 3^: 



(5.9) 



2fc5]Tre^"(*-*")(i?»(u,)(0jJ + 4(u,)(0jJ 



< ^^(1+^-0/2(1 + log l)(l-r)/2 I ^^^.g2a(t.-t„)||^jj|2_ 



1/2 



Proof. From (|5.2p . we observe that 

n 

2fc^Tfe2"(*'-*")i?'(u„)(0D 

i=l 

^ ^~'Y.( rf\^'^'^-'-\t-t,^^)\\uhu\\r-idty 

1 = 1 

Using p. Ill) , we find 

fc-'V( r r;/'e"(*'-*")(t-t,_i)l|u,t,||,_idt 



^ 1/2 r „ 



-I 1/2 



„2.(t.-.„)|,^.||2_ 



1/2 



(5.10) 



< 



,_i Jti_i L Jo 



'ti-i 



-,1/2 



[t~t,_^fdt 



It is now easy to calculate the remaining part for various values of r, s. For the sake of completeness, we present 
below the case when r = s = 0. 



^ / t^'^(t-ti^ifdt< / tdt + k"^^ t^'^ 



< iffc2(l + log-). 

k 

This completes the proof of the first half. For the remaining part, we observe from (|5.3p and (|4.3I) that 



(5.11) 2fc^- 



.s 2a(t.-t„) 



4(u,.)(0D < 



.2.(t.-*„)|,^.||2_ 



1/2 



4fcE (E / </^e"(*--'")(t-i,_i)/3(t. -t){5|lu,|l„+i + ||u,,|l,.+i}df)^ 
1=1 ,=i"'*j-i 



1/2 
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In Lemma 13.21 we find that the estimates of |luhtt||r-i and ||u/it||r+i are similar, in fact, the powers of ti are 
same. Therefore, the right-hand side of (|5.11|) involving ||uftt||r.-i-i can be estimated similarly as in (|5.10p . The 
terms involving ||u/i||t.+i are clearly easy to estimate. But for the sake of completeness, we provide the case, when 
r = s = 0. 

4'5'^E(E / e"(*--*")(t-t,_i)/?(i,-i)||Vu,|| dt) 
i=i j=i-^*j-i 

i=l j=l •'^j-i 

i=l •'0 Jo 



O 2 r2 ^ 



3g2(5fc 



This completes the rest of the proof. □ 

Lemma 5.2. Assume (A.\)-(A.2) and a spatial discretization scheme that satisfies conditions (Q\)-(R2) and 
fB4j. Let < a < min |(5, /iAi} , and 

which holds for < fc < fco, /cq > 0. Further, assume that u^(t) and V" satisfy anrf respectively. Then, 

there is a positive constant K such that 

n 

(5.12) UnW + e~2"*"fc^e2"*i|J?<iffc(l + log-), 

11 

(5.13) \\u\l + A;^{||^.ll^ + Iiat^jn<i^. 

Proof. For n — i, we put = in (j5.7p and with the observation 
we find that 

(5.14) dtlM' + Mml' + c^iim.i^) < -2i?Xu,0(^,)-24K)(|J. 
Multiply (|5.14p by ke^^^^ and sum over l<z<n<A^to obtain 

n— 1 n n 

\\ln\? - E(e'"' - 1)11^^11' + ^A'^E ll^^^ll' ^ -2fc5]e2"*'{ii;'(u,)($J +el(u,)(Cj} 

1=1 i— 1 z— 1 

n 

(5.15) < M E 11^^*11' + + log T)e'"*"+^ . 

Recall that v{t) = e"*w(t). Note that we have dropped the quadrature term on the left hand-side of (|5.14p after 
summation as it is non-negative. Finally, we have used Lemma |5. II for s = r = 0. We note that for < /c < fcg 

g2afe _ 

kXi 

and hence, 

2ak 1 ^ 1 

(5.16) ||i,f + ^-^^)kY, ||V|,f < Kk{l + log ^)e2"*"+^ 
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Multiply ()5.16|) by e to establish (|5.12p . Next, for n = i, we put </)^ = —Ah^^ in (|5.7p and follow as above to 
obtain the first part of (|5.13p . that is, 

n 

\\ui + kY,m\i<K. 
1=1 

Here, we have used (|5.9p for s = 0, r = 1 with a = replacing (f)]^ by A/i^j. 
Finally, for n — i, we put 0/j = 9t^j in (I5.7P to find that 

(5.17) 2m,f + Mlie.ll? < ~2aiql.{^),dt^,) - 2i?^K)(at|J - 24K)(9t|,). 

Multiply (|5.17l) by A:e^"*' and sum over 1 < i < n < A^. As has been done earlier, we can estimate the last two 
resulting terms on the right-hand side of (|5.17l) using (15.91) for r = s = as 

The only difference is that the resulting double sum (the term involving g*) is no longer non-negative and hence, 
we need to estimate it. Note that 

n n i 

(5.18) 2fc5]e2"*-a(g;(0,at|J = 27^^ ^ ^ e-(^-")(*-*^)a(|^., e"*'^*^,) 

i—1 i—1 j—1 



2 



i—1 i—1 J — 1 

Using change of variable and change of order of double sum, we obtain 

n i 2 

i=i j=i 

71 i i 

< if (7)fc E E e-(*-")(*-*^)) (fc E e-(^-")(*-*^) II A,,^^. | 

i=l j = l j = l 

<if(a,7)e(^-")'=fc2^fcE^"^'""^^*'"*'^ll^"^.7-|l'- 
i=i j=i 

Introduce I = i ~ j to find that 

n 

/<i<r(a,7)e('-"^'=fc'E^ E e-(^-")*M|A,i,_if for / - » - j 

4=1 l=i~l 



»i-i+l I 
i=l 1=1 

With change of summation, we now arrive at 

n n 

I < K{a, 7)e(^-")'=fc2 E ^ E e-^'-")*'- II A,i_,+i f 

n n — ^+1 

= if(a,7)e(''-")'=fc2^fc E e^(^-")*'-M|A,,|^.f for j = i - / + 1 
1=1 j=i 

n— 1 n 

(5.19) <X(a,7)e(^-")'=fc(fcEe"'*^"^*')(fcEll^''^jll') - ^• 



Combining ([5TTS)) - ([5TTg)) . wc find that 



n 7 ^ 

2fcEe'"*'a(g;($),5,$,) < ^ E 11^*^^ H' + ^■ 



2 

1=1 
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Therefore, we obtain 

n n—1 / 2ak i ^ 

(5.20) fc^e^"*'Ha,g,f + A^|||j|^<i^ + /ifc^ ^" ~ ' WIWI 

i=l i=l 

Use (15.121) and the fact that (e^"*^ — l)/fc < K{a) to complete the rest of the proof. □ 



Remark 5.1. We note that the restriction on k, that is < k < ko is not same in the Lemmas \4-l\ and \5.2\ 
Therefore, we take minimum of the ko 's from Lemmas \4-l\ awrf 1 5. 2\ and denote it as /cqO; then for all k satisfying 
< k < fcoO; 0,11 the result should hold. 

Analogous to the semi-discrete case, we resort to duality argument to obtain optimal L^(L^) estimate. Consider 
the following backward problem: For a given W„ and gi, let W;, n > i > 1 satisfy 

n 

(5.21) ((/.^,aiW,) -A*«(0^,W,) -i,)a(0/.,W,) = (0^, e^^t^g,), 0,, G J/.. 
The following a priori estimates are easy to derive. 

Lemma 5.3. Let the assumptions (A.2), (^Ij, (^2j and fB4) hold. Then, /or < k < k^, the following estimates 
hold under appropriate assumptions on W„ and g: 

n n 

\\Wo\\l + kY,e-'^'^{\\W,\\r+i + \\dtW,\\r-i}<K{\\WJl + kY,^^^^^ 

i=l j=l 

where r G {0, 1}. 

Lemma 5.4. Under the assumptions of Lemma \5.3\ the following estimate holds: 

n 

(5.22) e~'"*"fc^e2"*'l|$j2 < ^^2^ 

i=l 

Proof With 

W„ = (-A,)-i^„, g,=$, Vz 
we choose 4)^^ — in (|5.2ip and use (|5.7p to obtain 

n 

e'"'' lie.f = {^^, dtW,) - M^,,W,) - k ^ /3{tj - tMi^,^J) 

71 

= dt{^,,W,) - (9*^,, W,_i) - Ma(^., W,) - fc^/?(i, - i.)a(l„ W,) 

i 

= dt{^,,W,) + k{dtt„dtW,) + kJ2m - tj)a(|,, W,) + E\uh){W,) 

n 

(5.23) +4(u,)(W,)-fc5]/3(t,-i.)a(l.,W,). 

Multiply (j5.23p by k and sum over 1 < i < n. Observe that the resulting two double sums cancel out (change of 
order of double sum) . Therefore, we find that 

n n 

(5.24) kY^e'^'-Wa' + WU-i = fcE [fc(5t^.,5,W0 +i?'(u,)(W,) +4(u,)(W,)]. 

i=l i=l 

From (|5.2p . we observe that 

kY,E\nu){^,)<kY,i^ / (s-<.-i)||u,,.||_2||W,||2 

A 1 A 1 "J ti—-\ 



(5.25) <^e"'=(^*"e2"1|u,,,||V^)'^'(fcX:e~'"*-||W.||^)'^'. 
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Similar to (j5.16|) . we obtain 

(5.26) fc^e:.K)(W,)<if(fc3^ / e''^%\\n,,r + \M')ds) (fc J] g-^"*- ||W "^^ 

1=1 1=1 "^0 1=1 



and 



1/2 / JL s 1/2 



(5.27) kY,mt,,dtW,)<k{kY,e"'''\\dta') (A:^e-2"*i5tW,"2 

i—1 i—1 2—1 

Incorporating (|5.25p - (l5.27D in (I5.24p . and using Lemmas 13.21 and [Ol we find that 

n 

(5.28) kY,e'"'^m\' + WU-i < i^fcV"*". 

i=l 

□ 

Due to the non-smooth initial data, we need some intermediate results involving the "hat operator" which is 
defined as 

n 

(5.29) 0":=fc^0j, 

i=l 

This can be considered as discrete integral operator. We first observe, using (14.41) . that 

i i 



i — 1 St-^i St- ' 



j=i i=i 
Here 9^ means the difference formula with respect to i. Now rewrite the equations (j5.7p (for n = i) as follows: 

i 

(5.30) (9t|„ 0,J+A^a(|„ 0,J + dl[kY^p{U - t,)a(|,, = -E\un){ct^^,) - 

We multiply (|5.30p by k and sum over 1 to n. Using the fact that 9t^„ — observe that 

n 

(5.31) (9tt , 0J + Ma(i., 0J + 0J = {E\nn){ct,^) + 4(u„)(<AJ). 

i=i 

Lemma 5.5. Under the assumptions of Lemma \5/A the following estimate holds: 

n 

(5.32) |||„f + e-2"*"A:Ee'"'1|V|,|P<iffc^(l + log-). 

i=l 

Proof. Choose 4>j^ = in (I5.3ip for n = i, multiply by fce^"*' and then sum over 1 < i < n. We drop the third 
term on the left hand-side of the resulting inequality due to non-negativity. 

n n i 

(5.33) e^'^'-UA^ + ^,kY^e'-n'^U^ <kY.e'-''kY^{\E'{^h)^^^^^ 

i—1 i—1 i — 1 



From (|5.2p . we find that 



kj2w{^h){L)\<\{iz r (5-i,-i)iiu,,,ii_ids)iivi 
,=1 ^ .,=1^*^-1 
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Similar to the proof of Lemma 15.11 , we split the sum in j = 1 and the rest to obtain 
(5.34) 



kJ2 \EHu,,m)\ < Kk{l + ilog he-'^'^WVl 



n i n ^ 



i=l 



Therefore, 
(5.35) 
Similarly 

n i n ^ 

(5.36) k J2 e'"''fc E l^'aK)(|.)| < ^fc ^ e^"*' || V^.f + Kk'{l + log -)e' 

i=l — 1 i—l 

Incorporate (|5.35p - (|5.36p in (|5.33p to complete the rest of the proof. 
We are now in a position to estimate L°°(L'^)-norm of 
Theorem 5.1. Under the assumptions of Lemma lS.Sl the following holds: 



(5.37) 



" 1 
tn\\U\' + e-''''"kY,a,\\Vt,r<Ke{l + \og- 



where Oi — i^e^"*'. 

Proof. From (|5.7p with n = i and = (7^^^, wc obtain 

9*(^.||^.f ) - e2"^-{lll,:-if + ( 1 + 2A.a,||V|,;f 

(5.38) +2aM{i), I.) < -2E\xih){a,i,) - 2el(u^)(a,|,). 

We multiply (j5.38p by k and sum it over 1 < z < n to find that 



2cx.k 1 ^ 

f7„||C„f + (2/i - ^^^)fc^a,||Ve 



n-l 



(5.39) 



-2fc^a,a(g:,(0,l,)-2A:^£;'(u;,)(a.^,)-2fc^4(u;,)(a.^,). 



As earlier, using (|4.4p . we note that 

n n 71 z — 1 

(5.40) 2k aMQliO,^^) = 2A: ^ 7a(i, - 2fc ^ fc ^ 9t/3(i, - i,)a($^., a.^J. 

i=l 1=1 i=2 j=l 

The first term can be handled as follows (for some e > 0): 

n n n 

(5.41) 2k ^ ja{i,,a,^,) < efc ^ a, || V| + if (e, 7)^ ^ e^"*' || f . 

2—1 2—1 i — l 

For the second term, using similar technique as in (|5.19p . we observe that 

n i—l n 



(5.42) 2fc^fc^9t/3(t,-i,>(|^,a,^,) <efc^a,||Ve,f 

n i — l Sk—1 2 



i=2 j=l 
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2at„ 



Combining ([O0l) - ([O2l) . we find that 

n n n 

(5.43) 2k J2 <^^a{ql.{0, < skj^ V^.f + Kk^, e^"*' || V^, f . 

i—1 i—1 i—1 

From Lemma |5.1[ we obtain for r = and s = 1 

n n ^ 

(5.44) 2k J2 {E\uh){a,^,) + el(u;,)(r7,^ J} < efc^ V|,f + Kk^il + log -)e 

i=l 1=1 

Incorporate the estimates (|5.43p - (|5.44p in (|5.39p and choose e = to conclude 

2oLk 1 ^ ^ 

^n|l^JP + (M-^^^)fcE^^II^^'^ll'^^^'(l+l0gfc)e'"'"+^'^Ee'"*'ll^^«ll'- 
^ i=l i=l 

We multiply by e^"** and use Lemma 15.51 to complete the rest of the proof. □ 
We now obtain estimates of t] below. Hence forward, Kt means KTe^^ . 

Lemma 5.6. Assume (Al^, f-A-2j and a spatial discretization scheme that satisfies conditions (Bl), (B2^ and 
(B4.). Further, assume that U" and V" satisfy ^.6^ and i5.5]) . respectively. Then, for some positive constant K, 
there holds 

n 

(5.45) WvJ' + e-2"*"fc^ e^"*" ||r,,f < if,„fc(l + log -), 

n 

(5.46) I|r7j|? + e-^"*"fc^e2"*-||r7,||?<ift„. 

1=1 

Proof. We shall only prove the first estimate as the second one will follow similarly. For n — i, we put (pi^ = r/^ in 
(|5.8p . multiply by fce^"*' and sum over l<i<n<Nto obtain as in (|5.15p 

" (p2ak _ 1 \ " 

(5.47) ||f7„f + 2^^kJ2\\^'h^r < k ^ ^ ' IM' + 2A: ^ e^"*' A^^J. 

2=1 i=l i=l 



We recall from (|5.4p that 

(5.48) Alirj,) = -fe«,|„r;,J -fe(e„<,r;,) -fo(e„^„r7,). 
Using p.Sp and Lemma [521 we obtain the following estimates: 

&(e.,e„»7.) < lll.ir/'l|V|,f/^||Vr,,|| + ||V|,||||r,,||||Vr,,I| 

(5.49) < ell V,7.f + K\\rj^r + Kk^'^H + log ^)^/'|| V^,f 

(5.50) ^'(uLe.r?,) < £||Vr,j2 + X||V^J2 

(5.51) 6(e.,uLr,,) < £|| Vr,, f + X(|| V| + ||r,, f ) . 

Incorporate (|5.49p - ()5.5ip in ()5.48p and then in ()5.47p . Choose e = and once again use Lemma [5T2l Finally, use 
discrete Gronwall's Lemma to complete the rest of the proof. □ 



Remark 5.2. Combining Lemmas 1 5. 2\ and 1 5. (A we note that 

n 

(5.52) ||e„f + e-2"*"fc^ e^"*' ||ej2 < KtM.l + log 

i—1 

n 

(5.53) l|e„||? + e-2"*"fc^e2at.[jg^||2 

1=1 

Therefore, we obtain suboptimal order of convergence for ||e„||. 
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Below, we shall prove optimal estimate of ||e„|| with help of a series of Lemmas. 
Lemma 5.7. Under the assumptions of Lemma \5.6l the following holds: 



(5.54) \\rjJ\ + e-''^'-kJ2e'"''Hr<Kt„k^ 

i=l 

Proof. Put 4>f^ — e^"*'(— A;i)^^77j in (|5.8p for n = i. Multiply the equation by fce^"*'' and sum over 1 < i < n < iV 
to arrive at 

n n — 1 n 

(5.55) ||r)„f_, + II^'^II' ^ E(^^'"' - 1)11^^11-1 + 2fc^e^"*'AUe^"*'(-A,)-i,7.)- 

2—1 2 — 1 2—1 

From (15. 4p , we find that 

|2AU(-A,)-ir,J| < |26(e„<,(-A,0-ir,J 

(5.56) +biule,,{~Ah)-'v^)-bie^■.f^^A-^h)-'v^)\■ 
For the first term on the right hand-side of (|5.56|) . we use p.Sp to find that 

(5.57) |26(e„<,-A-i,7.)l<^l|edlll<l|i||77ji/'ll»7Jl'/'. 
Also, 

|26«,e„ -A-ir,,)| < |« . Ve„ -A-i^yJI + |« . V(-A-ir,,J,e,)| 

(5.58) < |« . Ve,, -A^'v.)\ + X||<||i||77j-'''ll^jr/'l|e.||. 

For Di — we note that 



2 r 

« • Ve„ -A-I77,) = ^ / <^,A(e.,,)(-A-i)r,,,^.dx 

2 

= / {A«,)e,,,(-A-i),7.,,+<,,e,,A((-A-i)r7,^^.)}rf2 



1.3 = 1 

(5.59) <X||<||i||e,||l!,7.l|i^i'l|r?.ir/'. 



Finally, from p.Sp . we find that 

(5.60) |2Ke.,e,,-A-V:)| <if||ed|(||ed|i + ||e,,||V2||e,||V2)||^j|V2||^^^ 

Now, combine (|5.56l) - (|5.60l) and use the fact that 

||e.||i< ||<||i + ||U'i|i<if 

to observe that 

|2AU(-A,)-ir,,)|<i^||e,||||,7j-/'l|r7jr/' 

(5.61) <i^liej|||r/,;||L/i'||r7,ir/2 + /f||»7.ll-/'ll»7.f/'. 
Incorporate (|5.6ip in (|5.55p and use kickback argument to obtain 



n n n 

2 



(5.62) ||r)„f_, + ^fc^ ||^,f < KkY,\mU + KkY^Wl 



i—l 2—1 2—1 



Finally, use Lemma 15.41 apply discrete Gronwall's lemma and multiply the resulting estimate by e 2"*^ to complete 
the rest of the proof. □ 
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Remark 5.3. From Lemmas \ 5.4\ and \ 5. 7\ we have the following estimate 

n 

(5.63) e-2"*"fc^e2"*'||e,f < Kt^ . 



i=l 



We need another estimate of rj similar to the one in Lemma 15.51 and the proof will follow in a similar line. For 
that purpose, we multiply (j5.8p by k and sum over 1 to n and similar to (j5.31l) . we obtain 

n 

(5.64) (9tr)„, 0J + fia{f,,^, </.,J + a{q'^{fj),<f>f^) = fc ^ 

1=1 

Lemma 5.8. Under the assumptions of Lemma \5.6\ the following holds: 

n 

(5.65) ||r,„ f + e-^"*" k J] e^"*' || Vr), f < if*„ k'{l + log -). 

1=1 

Proof. Choose — f]^ in ()5.64p for n = i, multiply by fce^"*' and then sum over 1 < i < n to observe as in ()5.33p 

n n i 

(5.66) e^"*" lirjjp + /xA: J] e^"*' || Vr)^^ < k J] e^-^'fc J] |AUj7.)l- 

2—1 i—1 j — 1 

We observe that 

i i 

(5.67) /c^|A^^(r,J| = fc ^ [&(<, e„ j),) + 6(e„ u^^, r,,) + 6(e,, e„ j),) 



Use dMl), (|53^ and ([535]) to obtain 



fc^e^"*'fc^|6(e„e„f,Jl</^A:^e2"*'fc^||e,||V2||e^.||3/2||v^^ 

2—1 2—1 i—1 J — 1 



i=l j = l j = l 

(5.68) <^t„fc'(l + log + £^ E e'"*' H^^^H"- 



Similarly, 

n i n '^1/9 

2=1 i—1 i—1 

n 

(5.69) <Xt„fc2 +£fc^ e^"*" II Vr),f. 

1=1 

and 



n i n i ^ 

kJ2e'^''kJ2 IKuLe.,0,)| < fcEe'"*'fcE (3 l((V • u?Je,>J| + !((< • V)0„e,)|) 

i^l j^l i^l j^l 

(5.70) <XA:^e2"*-(A:El|e.f) || V77. || < i^t„fc^ + efc E e'"*' II Vrj.f . 

i=l j=l i=l 

Combining these estimates, namely; (I5.68p - (l5.70p and putting e — /1/6, we conclude the rest of the proof. □ 
We present below a Lemma with optimal estimate for 77^. 
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Lemma 5.9. Under the assumptions of Lemma \5.6\ the following holds: 

n 

(5.71) tjrjj' + e-2"*"fc^ a,||r,j2 < Kt^k^l + log -). 

1=1 

Proof. We choose 0^ — ait]^ in (I5.8P for n = i. Multiply the resulting equation by k and sum it over 1 < i < n to 
find that 

n n — 1 n 

a„ ||r,„ f + 2/.fc ^ a, II f < (a)fc ^ 11^7, f - 2fc 5] a((;; (r,) , a,J7 J 

n 

(5.72) +2kY,K{<^^V^)■ 

i=l 

As in (I5.43p . we obtain 

n n n 

(5.73) 2fc^a,a(g;(,7),»7.)<efc^'T,||V,7.f 

1 i— 1 i— 1 

Following the proof technique leading to the estimate (|5.48p . we observe that 

n 71 71 

(5.74) 2fc^ AUa^J < efc^ a,;||Vr,,;f + /^fc ^ a,(|| V^,f + ||77,;f ). 

i—l i—1 i—1 

Substitute (15. 73^ - (15. 741) in (I5.72p and this completes the rest of the proof. □ 
Theorem 5.2. Under the assumptions of Lemma \5.b\ following holds: 

(5.75) ||e„|| <X^t-V2fc(i + log 1)1/2^ 



Proof. Combine the Lemmas 15.11 and 15.91 to complete the rest of the proof. □ 

Remark 5.4. We need not split the error e in ^ and r] in order to obtain optimal error estimate |5. 75[ ). However 
for optimal error estimate in L^-norm which is uniform in time, we need to split the error e„ = *7„ — Cn- 



6 Uniform Error Estimate 

In this section, we prove the estimate (j5.75p to be uniform under the uniqueness condition /i — 2iVt^^^ ||f ||oo > 0, 
where N is given as in p.l7p . We observe that the estimate (|5.37p involving ^„ is uniform in nature. Hence, we 
are left to deal with estimate of r7„. 

Lemma 6.1. Let the assumptions of Lemma \5.6\ hold. Under the uniqueness condition fi ~ 2iVt^^^ ||f ||oo > and 
under the assumption 

l + (^)fc>e^"^ 

which holds for < fc < fco, kg > 0, the following uniform estimate hold: 

1. 

k' 



3.1) ||»7„ II <i^r-i/2fc(l + log 1)1/2, 



where t„ ~ min{l,t„}. 

Proof. Choose 4>|^ = ij^ in (|5.8p for n = i to obtain 

(6.2) 5t||T7j' + 2/x||Vr,J+2a(g;(r,),r,J <2AUr7,). 

From ()5.4p . we find that 
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From the definition of N (see (|3.17p ). we note that 

(6.4) |6(r,„<,i7.)l <A^I|Vr,,f ||V<|1. 
Again with the help of (|3.8I) , we obtain 

|6(^,,<,,7JI + IKuLe.,^JI < ^ll^lll|V»7j||V<||i/2||A,<||i/2 

(6.5) <KTr'/\{l + loglnVri,\\. 

Since ||e,||2 < l|u*^l|2 + W^'h < Kt^^^'^, we conclude that 

(6.6) |&(e„e,,r7,)| < AVr3/4fc(i + iog^)i/2||Vr,J|. 
Therefore, from (|6.4p - (|6.6p . we find that 

(6.7) |AU»7JI <A^I|Vr,jn|V<||+ifrr3/4fc(l + log^)'/'||Vr,J. 



We recall from [23 that 



limsup ||Vuft(t)|| < ^l|f||oo, 



and therefore, for large enough z E N, say i > io we obtain from 

(6.8) |2AUe,)|<2iV^-i||f|U||V77.JP + ifrr3/'Ml + log^)^/'||V,7j. 
With ai ~ T^e^'^*S multiply f|6.2p by kai and sum over + 1 to n to obtain 

n n 

k J2 e'''''{^t\\v^\\' + M^V^f} + 2k J2 e^"*'a(g;(,7),»7.) 

2— zo + 1 2— zn + 1 

n 

(6.9) <2fc ^ <y^K{V^)■ 

i=io + i 

Without loss of generality, we can assume that iq is big enough, so that, by definition ti = 1 for i > io- We rewrite 
(HH) as follows: 

n n 

k e2"*Hat||r?,f + 2A*||Vr,,f} + 2fc^e^"*"a(g;(r,),r,,) 

i—iQ-\-l i—1 
n io 

(6.10) <^k J2 <^^Kiv^) + '2kY,^'"''<lliri\ri,)- 

We observe that the last term on the left hand-side of (IG.lOp is non-negative and hence is dropped. 

n— 1 n 

e'-'-WvJ'- E e"^'^{e"^''-l)\\rj,r+f,k ^ e'^^lVv^ 

i=io + l i=io + l 

+k E (M-27Vz.-i||f|U)e^-||Vr,j2<e2"*.o||r,,J|2 + 2fcEe2"*-g;(||Vr,,|l)||Vr,,|| 

i—iQ-\-l i—l 

(6.11) +Kk' J2 -■/'e^--(l + logi)3/^||V^, 
Under the assumption 



k' 

i=io + l 



l + (^)fc>e2"S 
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which holds for < A; < fco, ko > with < a < min {(5, we find that 

n n — 1 

|fc 5: e^-'llVr,,;!^- ^ e^-'(e2"'=-l)||r,,f 

z— io + 1 i— io + l 

(6.12) =fc E (f-^ll^)-^l|Vr7.f >0. 

i=io + l 

Due to uniqueness condition, we arrive at the following: 

n 

(6.13) k J2 (M-2A^i^-il|fl|oo)e2™||V77,||2>0. 

i=io + l 

Following the proof techniques of (|5.18p - (|5.19p . we obtain 

io io 

(6.14) 2A^e2"*^g;(||Vr,J)||V,7,:ll < i^fc ^ e^"*' |l V»7, f • 

1=1 i=l 

And 

n - n 

^ ,VV--(l + logi)V2l|Vr,J|<^fc 5: a.llVr,, 

i=io + l i=io + l 



n 

(6.15) +i^fc2(l+log-)A: V e2"*-rri/2. 
Incorporate (|6.12p - (j6.15p in (j6.1ip . use Lemma [^771 and (15. 63^ to observe that 

(6.16) +iffcf^e2"*-|lVr,J|2. 
Multiply by e~^"'' and under the assumption that 

(6.17) k e'"*- II V»7. f < ^toC^'(l + log r)- 



i=l 

we conclude that 

||r,J|<XC/2fc(i + log 1)1/2^ 

since zo > is fixed. Combining this result with (|5.37|) we complete the rest of the proof. 
We are now left with the proof (|6.17p . 

Lemma 6.2. Under the assumption of Lemma \5. 61 the following holds 

k f2 e'"*' II V77.IP < K,^,ti'k^{l + log I). 

i—l 

Proof. In (|5.47p . we use 

Kir],) = -bh{ule„rj^) -bh{e„V\rj^) 

<^||Vr,,f + X||e,f(||A,<|| + ||A,U^||), 
along with Lemma 15.71 and Theorem 15.21 to arrive at 

Jo 1 '0 

(6.18) HiciP +MfcE ii^^^ii' ^ ^t-o^' +^*.oVo'fc'(i +iog^)^Ee'"''^r 

i=l i=\ 

This completes the rest of the proof. 
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7 Conclusion 



In this paper, we have discussed optimal error estimates for the backward Euler method employed to the Oldroyd 
model with non-smooth initial data, the is, Uq G Ji. We have proved both optimal and uniform error estimate for 
the velocity. Uniform estimate is proved under uniqueness condition. The error analysis for the non-smooth initial 
data tells us that we need a few more proof techniques than the smooth case and proofs are more involved. 
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